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Abstract. New solar models are presented, which have been computed with the 
most recent physical inputs (nuclear reaction rates, equation of state, opacities, mi- 
croscopic diffusion) . Rotation-induced mixing has been introduced in a way which in- 
cludes the feed-back effect of the /x-gradient induced by helium settling. A parametriza- 
tion of the tachocline region below the convective zone has also been added in the 
computations. The sound velocities have been computed in the models and compared 
to the seismic sun. Our best model is described in some detail. Besides the new 
physical inputs, the most important improvement concerns the computations of 
/i-gradients during the solar evolution and their influence in slowing down rotation- 
induced mixing. This process can explain why lithium is depleted in the present 
Sun while beryllium is not, and meanwhile why ^He has not increased at the solar 
surface for at least 3 Gyrs. 
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1. Introduction 

For a long time the effects of microscopic diffusion were neglected in the 
construction of solar models. During the last two decades helioseismol- 
ogy has strongly constrained our view of the Sun and the importance 
of microscopic diffusion has been demonstrated by many authors (e.g. 
Bahcall et al., 1995; Christensen-Dalsgaard et al., 1996; Richard et al., 
1996 (RVCD); Brun et al., 1998a; Turcotte et al., 1998; Schlattl, 2002; 
Schlattl and Salaris, 2003). Microscopic diffusion is now systematically 
introduced in standard models. 

However some significant differences do remain between the he- 
lioseismic and theoretical standard models in the region below the 
convective zone. Several solutions for this discrepancy have been tested, 
such as opacity changes (Tripathy and Christensen-Dalsgaard, 1998) or 
turbulent mixing (RVCD; Brun et al., 1998, 1999 and 2002). Including 
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mixing below the convective zone is also necessary to account for the 
lithium depletion observed at the solar surface. This mixing however 
must not extend too deep as beryllium is not depleted (Balaehandran 
and Bell, 1997) and the ^He value has not increased significantly during 
the last 3 Gyrs (Geiss and Glocker, 1998) In previous papers, rotation- 
induced mixing was introduced in the computations but the feed-back 
effect of /Lt-gradients on the mixing was neglected or extremely simpli- 
fied. In RVCD, all macroscopic motions were supposed to disappear 
for a given critical /Lt-gradient, which was adjusted to satisfy the chem- 
ical constraints. Recently, the effects of diffusion-induced /v-gradients 
on the meridional circulation and their consequences have been more 
deeply studied (Vauclair, 1999; Vauclair and Theado, 2003; Theado and 
Vauclair, 2003 a and b). Numerical tools have been built which can be 
introdTiced in the computations of stellar evolution and applied to the 
solar case. 

We have computed new models with the Toulouse-Geneva Evo- 
lutionary Code, including the most recent physics for the opacities, 
equation of state, nuclear reaction rates, etc. We have also introduced 
diffusion and rotation-induced mixing in a consistent way, including 
the effects of /U-gradients. 

In section 2, we discuss the observed solar parameters and their 
uncertainties; the new physical inputs are presented in section 3 ; sec- 
tion 4 is devoted to the constraints on the mixing processes involved 
(rotation-induced mixing and tachocline); the mixing theory is reviewed 
in section 5 and the results, with our best model, are presented in 
section 6. 



2. The observed parameters 

The solar luminosity Lq, radius i?Q, and age allow one to constrain 
the initial Helium abundance, Yq, and the mixing length parameter, a, 
of the standard model. The luminosity Lq is obtained from the solar 
constant measured by satellite and varies with the solar activity. We 
take 3.8515 x 10^^ ergs s~^ (Guenther et al., 1992), which is the mean 
value between the results of the two space missions ERS on Nimbus 
7 and ACRIM on SMM. The slightly greater value used by Bahcall & 
Pinsonneault (1992) and the smaller value used by Bahcall et al. (1995) 
lie inside the error bars (see Table 1). 

Some recent papers have shown that the solar radius value of (6.9599± 
0.0007) X 10^° cm (Allen, 1976) was overestimated. Two different tech- 
niques now lead to a smaller value for the radius. The first one uses 
solar f-modes: Antia (1998) and Schou et al. (1997) have respectively 
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obtained 6.9578 x 10^° cm and (6.9568 ± 0.0003) x 10^° cm. The sec- 
ond technique uses the observation of the angular diameter of the 
Sun and atmosphere models. Brown & Christensen-Dalsgaard (1998) 
have found a smaller radius of (6.95508 zt 0.00026) x lO^'^ cm. For 
the calibration of our best model we take the mean observed value 
(6.95749 ± 0.00241) x 10^° cm, which is very close to the helioseismic 
value obtained with the f-modes. 

The last parameter used in the calibration is the solar age. Few 
determinations of the solar age have been done. The meteoritic de- 
termination obtained by Guenther et al. (1992) and by Wasserburg 
(See Appendix of Bahcall et al., 1995) is respectively 4.52 zt 0.04 Gyr 
and 4.57 lb 0.02. Using the helioseismic frequencies Guenther & Demar- 
que (1997) find 4.5 ± 0.1 Gyr and Dziembowski et al. (1999) obtain 
4.66 ±0.11 Gyr. 

The precise radius at the bottom of the convective zone is also given 
by the seismic sound speed profile. Christensen-Dalsgaard et al. (1991) 
found 0.713 lb 0.003. This value was confirmed by Basu & Antia (1997) 
with a smaller uncertainty; they gave : 0.713 zb 0.001, when they used 
only their models with a smooth composition profile at the bottom of 
the convective zone. 

The helioseismic inversion can also give the helium abundance in the 
convective zone (Dziembowski et al., 1991; Kosovichev, 1995; Basu and 
Antia, 1997; Richard et al., 1998). The most important uncertainty on 
the determination comes from the equation of state. With the OPAL 
equation of state (Rogers et al., 1996) the typical value obtained for 
the helium mass fraction Y is 0.249 ±0.002; with the MHD equation of 
state (Dappen et al., 1998; Hummer and Mihalas, 1988; Mihalas et al., 
1988) the results differ more from author to author (0.246 for Basu and 
Antia, 1997, 0.242 for Richard et al., 1998) due to the different inversion 
method. Basu et al. (1999) have recently studied the equation of state 
on the upper convective zone. They have found better agreement for the 
MHD equation of state than for the OPAL one on the helium ionization 
zone. However more detailed studies need to be done to reduce these 
uncertainties. 

The photospheric abundances of heavier elements used as constraints 
on our solar models are those determined by Grevesse &: Noels (1993). 
The values used to calibrate our models are summarized in Table L 



3. New physical inputs 

The present models have been computed including the most recent 
improvements in the physical inputs : 
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— Equation of state : Considering the high precision presently achieved 
in hchoseismology, it becomes necessary to take into account the 
relativistic effect on the electrons in the core of the solar model 
(Elliot and Kosovichev, 1998). In the present computations, we 
have introduced the new OPAL equation of state (Rogers and 
Nayfonov, 2002) which takes this effect into account. 

— Opacities : We use the latest OPAL opacities (Iglesias & Rogers, 
1996) completed by the Alexander &; Ferguson (1994) low temper- 
ature opacities. 

— Nuclear reaction rates : We have updated the nuclear routines 
using the NACRE compilation of nuclear reaction rates (Angulo 
et al., 1999) and added Bahcall's screening routine. 

— Microscopic diffusion : we have improved the computation of the 
microscopic diffusion velocity by using equation 4.7 of Montmerle 
& Michaud (1976) for all elements other than helium. Turcotte et 
al. (1998) have pointed out a difference in the surface abundance 
variations of C, O, and other metals (Z), during the evolution, 
between the Toulouse and Montreal solar models (see table 3 of 
Turcotte et al., 1998). This difference was mainly due to the ther- 
mal part of the diffusion velocity which was underestimated in the 
expression used by RVCD. With the Montmerle &i Michaud ex- 
pression, the Toulouse models have very similar surface abundance 
variations of C, O, and Z as those of the Montreal models 



4. Constraints on mixing processes 

4.1. Chemical constraints on mixing processes 

Models including helium and metal diffusion cannot reproduce the 
observed lithium depletion (1/140). Although part of the lithium deple- 
tion may occur during the pre main-sequence phase, the main-sequence 
lithium depletion is likely to be very much larger than that induced by 
microscopic diffusion alone. Lithium diffuses at about the same rate 
as helium, so that its abundance would decrease by some 10% only 
by diffusion. On the other hand, mixing can bring matter from the 
outer layers down to the place where lithium is destroyed by nuclear 
reactions, which is much more efficient. Moreover diffusion builds a 
steep helium gradient below the convective zone which leads to a spike 
in the sound velocity compared to the seismic sun in this region. These 
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Figure 1. Location of lithium and beryllium nuclearly-depleted region with time 
(Myr) in a standard model including diffusion. 

results suggest the presence of some mild mixing below the convective 
zone of the internal Sun. 

Macroscopic motions are needed to bring up to the convection zone 
the nuclearly lithium-depleted matter and to smooth the helium-gradient 
at the top of the radiative interior. On the other hand we now know that 
beryllium is probably not depleted in the Sun and that the ^He/^He 
ratio cannot have changed by more than 10% during these last 3 Gyrs 
(Geiss and Gloeckler, 1998). This means that the mixing must not 
bring too much helium to the convective region and that the beryl- 
lium nuclearly-depleted region must not be connected to the convective 
region through the mixing. 

This put strong constraints on the mixing : it must be deep enough 
to deplete lithium and to smooth the helium-gradient below the con- 
vective region but it must be also shallow enough to keep constant 
the beryllium abundance and the ^He/^He ratio during the last 3 
Gyrs. We have represented on figure 1, the lithium and beryllium 
nuclearly-depleted regions as a function of the age in a standard model 
(including diff'usion). According to the location of the lithium and 
beryllium nuclearly-depleted regions, the mixed zone may extend on a 
shallow region below the convective zone, down to 0.6i? at the beginning 
of the solar life (where R represents the radius of the sun at that time) 
and down to O.55i?0 at the solar age. 
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4.2. Observational constraints on the tachocline 



Inversions of rotational splittings have shown that a strong horizontal 
shear exists in the solar rotation profile at the base of the convection 
zone in the so-called solar tachocline. Several attempts at estimating 
the characteristic thickness of this region have been made from helio- 
seismic inversions. Values of this thickness range from 0.04 — O.O5i20 
(Basu, 1997; Antia et al., 1998; Basu and Antia, 2001; Charbonneau 
et al., 1998, 1999; Corbard et al., 1998) to O.O9i?0 (Kosovichev, 1996; 
Kosovichev et al., 1998), the thickness being defined as the region where 
the rotation rate variation inside the tachocline increases from 0.08 
to 0.92 of its maximum value. More recently Corbard et al. (1999) 
concluded that the tachocline width is very likely below O.O5i?0 and 
Elliot k Gough (1999) obtained a value of O.O19i?0. 

As previously mentioned, the beryllium and ^He/^He solar obser- 
vations favor of a shallow mixed region below the convective region. 
Therefore the tachocline may provide an adequate mixing. However 
from helioseismic observations the tachocline is extremely thin and, 
supposing that its width has not been strongly increased in the past, 
the tachocline-induced mixing is not deep enough to produce a smooth 
helium gradient below the convective zone. Brun et al. (1999) have 
then considered a variable tachocline width during the solar evolution. 
In these conditions they have computed models in agreement with 
helioseismic and chemical observations. But there is no evidence and 
no explanation for a variable tachocline with time. 

These results are in favor of extra mixing, probably related to the 
rotation-induced thermal imbalance. In the following we present solar 
models including both the effects of a tachocline and a rotation- induced 
mixing. In section 5 we present the non-standard physical inputs in- 
troduced in the computations : we review the rotation-induced mixing 
theory together with the main equations, we discuss the stabilizing 
effects of /x-gradients and we present the chosen tachocline modeling. 
Section 6 is devoted to the computational results. 



5. Mixing theory 
5.1. Rotation-induced mixing 

The theory of rotation-induced mixing in the presence of //-gradients 
is introduced in the models. 

When the rotation depends little on latitude, the vertical component 
of the meridional velocity can be expressed in term of the Legendre 
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function of order 2 : 

u{r,e) =Ur{r)P2{cos9) (1) 

r being the radial coordinate and 9 the colatitude. The horizontal com- 
ponent of the meridional velocity also depends on the vertical velocity 
amplitude Ur{r) : 

ug = ^ ( pr^ Ur) sin 6 COS (2) 

2pr dr 

In the case of the stationary regime and negligible differential rota- 
tion (as observed from helioseismic studies) the vertical amplitude can 
be written : 

P L 

E^, gathers the terms related to the /U-gradients while the classical terms 
appears in Eq : 



E = ^ {- — 



{ r d \ dk , , . 1 Ht .1 



Here p represents the density average on the level surface (~ p) 
while Pm is the mean density inside the sphere of radius r; Ht is the 

temperature scale height; A represents the horizontal ji fluctuations ^; 

Xn and xt represent the derivatives: 

Traditionally the meridional currents induced by the Eq_ term arc 
called r^-currents while the ones induced by are called /x-currents. 
In the most general case, ^-currents are opposite to J7-currents. 

The rotational mixing is then introduced in the computations follow- 
ing Zahn(1992) and Maeder k, Zahn (1998). The meridional circulation 
is supposed to induce shear flow instabilities which, in a density strati- 
fied medium may lead to large horizontal diffusivities. The combination 
of the meridional circulation and horizontal mixing may be treated as a 
vertical effective diffusion process described with an effective diffusion 
coefficient. The chemical transport for the mean concentration of a 
species i is then expressed : 



rtv04.tex; 2/02/2008; 8:58; p. 7 



8 



Richard, Theado and Vauclair 



where 

J. _ „ , [rUr{r)f 

^turb — -L>v -\ 5777:; [o) 

6UL>h 

Dy and Dh represent respectively the vertical and horizontal part of the 
shear-induced anisotropic turbulence due to the transport of angular 
momentum. 

In the stationary case, according to Chaboyer & Zahn (1992) : 

In the present paper, the horizontal diffusion coefficient is approxi- 
mated following Maeder & Zahn (1998) : 

Dh = Chr\Ur\ (10) 

The efficiency of depends on the Richarson criterion (Maeder, 1995; 
Talon and Zahn, 1997). Here we chose to parametrize this coefficient 
in terms of (rUr) as for the vertical diffusion coefficient : 

= Cyr\Ur\ (11) 

In this case, A can be written : 

A=- — (12) 
6Ch\Ur\ dr ^ ' 

and we can also write : 

Dturb = (^Cy + -^^Q^^Prl = aturbr\Ur\ (13) 

with aturb = + ^ . Cy and Ch are then unknown parameters. 

According to the assumption of strong anisotropic turbulence they must 
satisfied the condition : 

Cy « Ch (14) 



5.2. Stabilizing effects of /x-gradients 

It has been known for a long time that mixing processes in stars may 
be stabilized in the regions where the mean molecular weight rapidly 
decreases with increasing radius (Mestel, 1965; Huppert and Spiegel, 
1977). This occurs specifically in the nuclear burning core. The im- 
portance of the feed-back effect due to nuclearly-induced /x-gradients 
was well recognized (see Mestel and Moss, 1986), however the feed- 
back effect due to the diffusion-induced /U-gradients on the meridional 
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circulation was not included until recently. In previous work on solar 
models, it was inferred that the rotation induced mixing becomes in- 
efficient as soon as the /v,-gradient becomes larger than some critical 
value (RVCD). This value was adjusted to obtain the lithium depletion 
observed and to satisfy the constraint on the ^He/^He ratio. 

Here we include a more precise treatment of the physics of diffusion- 
induced /U-gradients and the diffusion-circulation coupling as discussed 
by Vauclair & Theado (2003) (hereafter referred as paper I), Theado &; 
Vauclair (2003a, paper II) and (2003b, paper III). While the /^-terms 
have definitely to be taken into account in the computations, the details 
of the coupling process which occurs when these term becomes of the 
same order of magnitude as the classical terms remain difficult to han- 
dle. In paper I, we introduced an analytical approach of the processes 
and discussed an approximate solution in a quasi-stationary case. We 
showed that the coupling between helium settling and rotation-induced 
mixing can slow down both the mixing and the settling when the ji- 
terms are of the same order of magnitude as the other terms. In paper 
II, we presented the results of a 2D numerical simulation of the con- 
sidered processes which helps visualizing the situation. In the example 
which has been chosen (slowly rotating low mass stars), the /x-terms 
become of the same order of magnitude as the classical terms in a short 
time scale compared to the main-sequence lifetime, thereby creating a 
frozen region where the circulation does not proceed anymore. How- 
ever when this occurs the equilibrium between the opposite meridional 
currents is permanently destabilized by the helium settling below the 
convective zone. As a consequence, a new motion develops, which mixes 
up the zone polluted by diffusion. Below this region the circulation 
remains nearly frozen. This process acts as a restoring system for the 
element abundance variations and increases significantly their settling 
time scales. We then suggested that diffusion and mixing react in such 
a way as to keep both the horizontal and vertical ^[/-gradients constant 
in the frozen region, while they proceed freely below. 

5.3. Tachocline 

Macroscopic motions due to the tachocline can be introduced in the 
computations by adding a diffusion coefficient Dfacho in the equation 
for the time evolution of the concentration of chemical species. In this 
case the turbulent diffusion coefficient -Djurfe which appears in equation 
7 must be replaced by an effective diffusion coefficient -De// '■ 

^eff = Dfurb + Dtacho (15) 
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Spiegel & Zahn (1992) and Brun et al. (1998b) have obtained an ex- 
ponential diffusion coefficient for Dfacho- Here the tachocline is treated 
as follows : 

Dtacho = Dcz exp(ln 2 — (16) 

where D^z and r^z are respectively the value of Dtacho a-^d the value of 
the radius at the base of the convective zone. A is called the half width 
of the tachocline. D^z and A are free parameters. Here the absolute 
size of the tachocline is assumed constant with time. 



6. Computational results 

The diffusion-circulation coupling and the tachocline described in sec- 
tion 5 have been introduced in the computations of solar models. 

6.1. Constraining the free parameters 
6.1.1. Rotation-induced mixing 

The free parameters of the rotation-induced mixing are Ch and Cv 
(or aturb) which determine the efficiency of the turbulent motions. The 
competition between mixing and diffusion then determine the construc- 
tion of /Lt-gradicnts and therefore also the strength of ^-currents. 

In particular the //-currents strongly depend on Ch which is directly 
related to the horizontal turbulence : a strong horizontal turbulence 
tends to homogenize the horizontal layers and thereby smooths the 
horizontal jU-gradients. In this case the small induced jU-currents may 
not be able to compensate for the 0-currents and the mixing may 
remain very efficient during all the solar life. On the contrary, weak 
horizontal turbulence leads rapidly to important jU-currents which may 
strongly reduce the mixing. 

In the computations, we try to adjust the mixing free parameters to 
produce a mixing both : 

— efficient and deep enough to smooth the diffusion-induced helium 
gradient below the convective zone, 

— weak and shallow enough to prevent beryllium from nuclear de- 
struction in the external layers of the Sun. 
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6.1.2. Tachocline 

By analogy with the usual definition of the tachocline (of section 5.3) 
we define the tachocline thickness uj as the region where the effective 
diffusion coefficient Dtacho (cf eq. 16) increases from 0.08 to 0.92 of 
its maximum value Dcz- As a result the constraint on u deduced from 
helioseismic observations (w < O.O4i?0) leads to a constraint on A 
which is : 



Considering this constraint, the two free parameters of the tachocline 
(A and Dcz) are adjusted to reproduce the observed solar lithium deple- 
tion without depleting beryllium. For that we will assume that lithium 
depletion during the pre-main-sequence is negligible. As a consequence, 
the deduced values for the tachocline parameters will be upper limits 
since lithium may be partly depleted on pre-main-sequence. 

6.2. Results 

The characteristics of our best model are given in table I. This model 
is calibrated to reproduce at the solar age the observed luminosity and 
radius of the Sun. This calibration is achieved by adjusting the two main 
free parameters of the code which are the mixing length parameter a 
and the initial helium abundance Yq. The obtained values are : 



As in RVCD, the rotation is assumed uniform in the models, but 
the rotational velocity decreases with time according to a Skumanich 
(1972) law, i.e. O oc 2 . The coefficient is adjusted for an initial surface 
velocity Vrq =100 km s~^ and a surface velocity at the age of the Sun of 
\rq=2 km s~^. In our computations, the rotation- induced mixing and 
the tachocline parameters have been chosen to reproduce the smooth 
diffusion-induced helium-gradient below the convective region and the 
observed lithium depletion while keeping the beryllium abundance con- 
stant. Their values are consistent with those proposed by Theado and 
Vauclair (2003b) for solar type stars. 

This model is computed with the following values of the rotation 
induced mixing parameters : 



A < 0.01135i?( 



'0 



(17) 



a = 1.753739 and Yq = 0.268039 



Ch = 7000 



and aturb = 1 



and the following tachocline : 



Drz = 2.5 X 10^ and A 



0.048 X 10^° = 6.900 x lO'^R, 
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Table I. Constraints on solar models and characteristics of our best 
model. 





Constraints 


Model 


Age (Gyrs) 


4.57±0.02 


4.606 


M (g) 


(1.9891±0.0004) 10^^ 


1.9890 X 10*^ 


L (erg s ) 


(3.8515±0.0055) 10^^ 


3.8523 X 10*^ 


R (cm) 


(6.95749±0.00241) 10^° 


6.95654 X 10^° 


rbzc/-R© 


0.713±0.003 


0.718 


F.uri (EOS OPAL) 


0.249±0.002 


0.2507 


(Z/X)0 


0.0245 


0.0245 


Li/Lio 


1/140 


1/133.8 


Be/Beo 


1 


1/1.219 


^Hc/*He 


^He/*He (ti t©) < 10% 




ti = 1.6 Gyrs 




1.614 X 10"" 


tg = 4.6 Gyrs 




1.642 X lO"'' 



Both the conditions Ch » Cy and A < O.OllSSi?© are satisfied. 

Figure 2 displays the hehum gradient of our model below the convec- 
tive region and it compares its sound speed square with the one deduced 
from helioseismic observations (Basu et al., 1997). For comparison the 
figure also shows the results obtained for a standard model includ- 
ing element segregation. Thanks to the mixing the helium-gradient is 
strongly smoothed and as a consequence the spike in the sound velocity 
observed in the case of pure element settling is reduced. However in 
the deep interior the rotation-induced mixing does not improve the 
agreement between models and helioseismology. 

Figure 3 shows the \Eq\ and \E^j\ profiles inside the star at differ- 
ent evolutionary stages. The sun arrives on the main-sequence with 
nearly homogeneous composition. At that time \E^\ is much smaller 
than \Eq\ everywhere inside the star. The meridional circulation and 
element settling can take place in the radiative region of the star. Just 
below the convective zone, under the effects of the tachocline material 
is strongly mixed up and homogenized with the convective envelope. As 
a result in this region, the /i-gradients cannot form and the //-currents 
remain negligible. On the other hand below the tachocline region, as a 
consequence of the helium settling and the f2-currents, the //-currents 
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Figure 2. Results at 4.606 Gyrs : the solid lines represent the results of our best 
model, the dashed lines show the results of a standard model including element seg- 
regation. Left panel ; He profile below the convective zone. Right panel : comparison 
between (sound speed square) deduced from helioseismology (Basu et al., 1997) 
and the computed one. 



increase, thereby reducing the mixing efficiency in a non-negligible way 
as described in the previous section. They also increase in the core 
because of nuclear burning but other physical processes can occur there 
which are not taken into account in the present computations (e.g. mild 
mixing in the internal layers at some stages of the solar evolution). 

Figure 4 displays the diffusion coefficients {D^ff, D^urb and Dtacho) 
inside the star at different evolutionary stages. The tachocline-induced 
mixing is significant over a narrow region located above O.6i?0 which 
means above the beryllium nuclear destruction region. Below the tachocline 
region, the mixing is due essentially to the rotation-induced mixing. In 
this region because of large //-currents, the diffusion coefficient becomes 
very small during the solar evolution. 

Finally at the beginning of the solar life, a strong mixing takes place 
in the whole radiative interior. Just below the convective zone, in the 
tachocline region, this mixing remains very strong during all the solar 
evolution. Below the tachocline region, the meridional circulation be- 
come rapidly very weak (after 2 Gyrs) . As a result the convective region 
is rapidly disconnected from the beryllium-nuclearly depleted region 
while it remains connected to the lithium one through the tachocline 
during all the evolution. 

Figure 5 displays the lithium and beryllium depletion with time at 
the solar surface. At the age of the Sun the lithium is depleted by a 
factor 133.8 while the beryllium depletion is lower than 18%. 
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Figure 3. Evolution of Q, and /j-currents with time inside the Sun. The graphs show 
the variations with depth of both \Eq\ and \E^\. 

Therefore the chosen values of the free parameters reproduce quite 
well the lithium depletion observed in the Sun and lead to a small 
destruction of beryllium. Table I also shows that the constraints on 
the ^He/^He ratio is satisfied as well as the constraints on the initial 
helium abundance ("Kj„r/=0.2507) which is consistent with the value 
deduced from hclioseismology by Basu & Antia (1997) and Richard et 
al. (1998) : Yq^^utJ = 0.249 ±0.002 (with the OPAL equation of state). 
The position of the convcctivc region is however slightly overestimated. 

The lithium abundance in the Sun may be compared to that ob- 
served in galactic clusters of known age. Table II then compares the 
lithium abundance in our model with that observed in four of these 
galactic clusters. The parametrization of the mixing and the tachocline 
below the convective zone correctly reproduces the lithium abundance 
evolution with time. 



7. Conclusion 

The chemical composition of the solar interior strongly relies on the 
competition between several physical processes : microscopic diffusion, 
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Figure 4- Evolution of the diffusion coefficients with time inside the Sun. The dot- 
ted linos show the turbulent diffusion coefficient Dturb, the dashed lines show the 
tachocline-induced diffusion coefficient Dtacho and the solid lines show the effective 
diffusion coefficient Def f=Dturb+Dtacho- 
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Figure 5. Lithium and beryllium abundance variations with time (yrs) at the solar 
surface. 
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Table II. Lithium abundances observed in 
galactic clusters and lithium abundances of our 
model 





Observations 


Model 


Age 


Cluster 


A(Li) 


A(Li) 


0.05 Gyrs 


a Persei 


3.2±0.3^ 


3.25 


0.1 Gyrs 


Pleiades 


2.9±0.3 ^ 


3.17 


0.8 Gyrs 


Hyades 


2.5±0.3 3 


2.55 


1.7 Gyrs 


NGG 752 


1.7±0.5* 


2.01 



^ Balachandran ct al. (1996) ^ Thornburn et al. (1993) 

^ Soderblom et al. (1993) " Balachandran (1995) 



rotation-induced mixing and tacliocline mixing. Up to now, the feed- 
back effect of helium setthng on the mixing was not taken into account 
in the computations, although the necessity for the mixing to stop 
at some depth below the convective zone appeared clearly from the 
observations of lithium, beryllium and helium 3. For this reason, RVCD 
had introduced an artificial cut-off of the turbulent diffusion coefficient 
for a given |Lt-gradient adjusted as a parameter. The present computa- 
tions are more precise as they take into account the coupling between 
microscopic diffusion and rotation-induced mixing as described in pre- 
vious papers by Theado and Vauclair 2003 a and b. Mixing is naturally 
slowed down and stopped when the so-called /x-currents arc opposite 
to the ri-currents. The present modeling procedure also includes the 
most recent physical parameters for the opacities, equation of state, 
and nuclear reaction rates. The seismic sun is nicely reproduced by 
our best model, although some discrepancy by about 0.5 percent still 
remains. The same modeling procedure applied to solar-type stars in 
galactic clusters also reproduces correctly the lithium depletion with 
time. 
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